Correlation effects in p-electron magnets: the case of Rb02 
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We present results of GGA+[/ calculations for the magnet" Rb02, where magnetic properties 
are due to partially filled oxygen p orbitals. We show that on-site interactions on the oxygen sites lead 
to a strong tendency towards the formation of an orbit ally polarized insulating state, in contrast 
to the half-metallic behavior predicted for this class of compounds within pure LDA/GGA. The 
obtained energy differences between different orbit ally ordered configurations are sizeable, indicating 
an orbital ordering temperature higher than the antiferromagnetic Neel temperature of ~15 K. Our 
results demonstrate the importance of correlation effects in p electron magnets such as Rb02. 



Recently, various cases of "d^-magnetism" have been 
reported where, in contrast to the more conventional case 
of partially filled d (or /) electronic shells, magnetic prop- 
erties arise from partially filled p orbitals [1]. In most 
cases, the corresponding phenomena are defect-induced, 
and systematic studies are hampered by poor repro- 
ducibility and a wide spread in experimental data. Fur- 
thermore, theoretical investigations based on first prin- 
ciples density functional theory are plagued by the well- 
known deficiencies of the local spin-density and general- 
ized gradient approximations (LDA and GGA) [2]. To 
further explore p electron magnetism as alternative op- 
tion for spintronic applications, it is therefore desirable to 
study intrinsic p electron magnetism in pure, i.e. mostly 
defect-free, bulk materials, which can provide important 
insights in the underlying mechanisms and can serve as 
benchmarks for currently used theoretical approaches. 

Here, we discuss the case of rubidium superoxide, 
Rb02, which is part of the family of alkali-superoxides 
AO2 {A = K, Rb, or Cs) 011]. The magnetic properties 
of these systems result from the partially filled p electron 
levels of the superoxide anion . As shown in Fig. [l^, 
the hybridization between atomic p orbitals within the 
O2 units gives rise to bonding (a^, TVx^y) and anti-bonding 
(cr*, TT*^) molecular orbitals (MOs). The nine p elec- 
trons per occupy all bonding states, but leave one 
hole in the antibonding tt* ^ orbitals (see Fig. [l^), which 
can give rise to spin polarization and magnetic order (see 
Fig.[T]3). Indeed, below the Neel temperature Tn ~ 15 K, 
Rb02 orders antiferromagnetically, with parallel orien- 
tation of the magnetic moments within the tetragonal 
(001) planes and antiparallel orientation between adja- 
cent planes [3l U . 

Rb02 also undergoes several crystallographic phase 
transitions [4, 5 , which involve various small distortions 
from an "average" high symmetry structure of the CaC2- 
type (space group /4/mmm, see Fig. [T]). This structure 
consists of two interpenetrating body-centered tetragonal 
lattices of Rb+ cations and anions with the molecular 
0—0 bond oriented parallel to the c axis. For simplicity, 
and since the full space group symmetries of the various 
low temperature phases have not been fully established, 
yet, all calculations presented in this work are performed 




FIG. 1: (color online). Left: Averaged tetragonal structure of 
Rb02. Oxygen and rubidium atoms are represented by large 
(red) and small (blue) spheres, respectively. In the following, 
the two O2 molecules within the tetragonal unit cell will be 
denoted as O^^ and O'^^ Thick solid (green) lines indicate 
the hoppings included in the TB model. Right: (a) Energy 
levels of O^ molecules formed from atomic p states, (b) Spin 
polarization splits the majority and minority 7ix,y states. The 
degenerate half- filled minority-spin n^^y orbitals give rise to 
the half-metallicity obtained within LDA/GGA calculations, 
(c) An additional splitting results in an orbitally polarized 
insulating state consistent with experiment. 

for the average CaC2-structure shown in Fig[l] 

In the following, we present results of first princi- 
ples density functional theory calculations using the 
Quantum-ESPRESSO code package [6 , employing the 
GGA exchange-correlation functional of Perdew, Burke, 
and Ernzerhof [7 and Vanderbilt ultrasoft pseudopoten- 
tials [8 (USPP). We analyze the effect of on-site Coulomb 
repulsion on the electronic structure of Rb02 using the 
GGA+/7 approach [9, 10 , and demonstrate that Rb02 
exhibits a strong tendency to form an orbitally polar- 
ized insulating state (see Fig. ^) , in contrast to the half- 
metallic character obtained within a simple LDA/GGA 
calculation. The orbital order occurs without imposing 
any symmetry- lowering of the crystal structure. Our re- 
sults demonstrate the importance of correlation effects in 
p electron magnets such as Rb02. 

We start by performing a full structural optimization of 
Rb02 within tetragonal lA/mmm symmetry. These re- 
laxations are done for antiparallel alignment of the mag- 
netic moments of 02^^ and 02^^ (see Fig. [l]). We use 
a plane wave cutoff energy of 30 Ry and alOxlOxG 
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FIG. 2: (color online). DOS of Rb02 for FM (a) and AFM (b) 
configuration. Black line represents total DOS; striped areas 
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FIG. 3: (color online), (a) Diagonal elements of the occupa- 
tion matrix as a function of the Hubbard U parameter, (b) 
Evolution of the total DOS (black line) and projected DOS 
(vertical and horizontal stripes for 0{px) and 0{py), respec- 
tively) for different U values. The solid vertical line indicates 
the Fermi energy. 



/c-point grid, which result in convergence of the total 
energy better than 1 meV. We obtain lattice parame- 
ters a = 4.20 A and c = 7.07 A, in very good agreement 
with experimental values of 4.22 A and 7.00 A, respec- 
tively ^j. The calculated 0—0 bond length is 1.36 A 
which agrees very well with the experimental estimate of 
1.34 A [n]. 

The calculated densities of states (DOS) for both fer- 
romagnetic (FM) and antiferromagnetic (AFM) config- 
urations are shown in Fig. [2] It can be seen that in 
the plotted energy range around the Fermi level the 
DOS has almost exclusive oxygen p character, and ex- 
hibits a peak structure closely resembling the molecu- 
lar energy scheme depicted in Fig. [l] Thus, the rela- 
tively weak interaction between neighboring molecules 
leads only to a small broadening of the corresponding 
energy levels, i.e. the formation of very narrow bands 
that essentially retain the character of the corresponding 
MOs. The DOS for the two different magnetic configu- 
rations show only minor differences. In both cases, the 
Fermi level bisects the (local) minority tt* ^ band which 
leads to metallicity (half-metallicity for the FM case), 
similar to LDA calculations for KO2 [l2] and GGA cal- 
culations for rubidium sesquioxide, Rb406 [13l[T4]. The 
latter material, which contains a mixture of superoxide 
and nonmagnetic peroxide anions, was predicted to be a 
half-metallic ferromagnet with an estimated Curie tem- 
perature of 302 K [131 HI]- However, experiments in- 
dicate that Rb406 is indeed a magnetically frustrated 
insulator exhibiting spin-glass-like behavior, consistent 
with more recent calculations using hybrid functionals 
and LDA+t/ [151 H^. It was also pointed out in Ref. [15] 
that even Rb02, which is experimentally well known to 
be an insulating antiferromagnet, is predicted to be half- 
metallic within LDA. This is consistent with our GGA 
calculations, which favor the FM over the AFM con- 
figurations by an energy difference of 6 meV per for- 
mula unit (f.u.), thus suggesting a general inadequacy of 
LDA/GGA for the treatment of partially filled molecular 



states. 

It is well known that LDA/GGA often fail to repro- 
duce the correct insulating ground state for "strongly 
correlated" materials such as transition metal oxides and 
many / electron systems (see e.g. Ref. [10 ). This failure 
can be corrected within the (LDA/GGA) +/7 method, by 
introducing a correction term accounting for strong on- 
site Coulomb interactions 0lQj. For the implementation 
used in this work the corresponding energy correction has 
the form [iTllIH]: 




Here, cr^ and m indicate the site, spin, and orbital 
character, and the sum runs over all atoms/orbitals to 
which the Hubbard correction is applied, n^^/ is the or- 
bital occupation matrix obtained by projecting the occu- 
pied valence functions IV^^^) band index) onto atomic- 
like states 10^): 

nw = E4>(^LI^™)<^-'I^L) • (2) 

Here, f^^ is the occupation of the corresponding wave- 
function. In the following we suppress the superscript 
Rcr where possible; if not stated otherwise, occupation 
numbers refer to the local minority-spin p channel on the 
oxygen sites. 

To assess the effect of on-site electron correlation for 
the present case of a partially filled p electron system we 
now perform a series of GGA+/7 calculations with vary- 
ing Hubbard U parameter for the p orbitals on the oxygen 
sites. Within this scheme, an insulating state can in prin- 
ciple be achieved via orbital polarization, i.e. a preferred 
occupation of one of the degenerate tt* orbitals driven by 
correlation effects (see Fig. [l]^). In the following we only 
consider the AFM configuration, and in order to allow 
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FIG. 4: (color online), (a) Total energy of orbitally polarized 
states relative to the non-polarized state as a function of U. 
(b) Energy of different orbitally ordered states relative to the 
{x,y) state. All energies correspond to the simple tetragonal 
unit cell. 




rxM rzRA zxRMA r 



3 4 5 6 7 
f/(eV) 



FIG. 5: (color online), (a) Comparison of the band struc- 
ture for TT* bands: GGA (thick/black line) and TB model 
(thin/red line). The two MLWF per O2 site used as basis for 
the TB model are also shown, (b) Band gap as a function of 
U. Symbols for GGA+[/ data indicate the same orbital order 
patterns as in Fig. [3] 



the system to converge to an orbitally polarized state, 
we initialize the elements of the orbital occupation ma- 
trix accordingly [19 . The resulting converged values of 
the occupation matrix elements as a function of U are 
summarized in Fig. [3^. Note that since the occupation 
matrix in the current GGA+/7 implementation refers to 
atomic Px^y states, which give rise to both the fully oc- 
cupied bonding and the partially occupied antibonding 
MOs, a fully occupied tt* MO and fully unoccupied tt* 
MO correspond to atomic occupation numbers rixx = 1.0 
and n^^ = 0.5, whereas for the corresponding unpolar- 
ized state Uxx = ^yy = 0.75. It can be seen from Fig. [3^ 
that no orbital polarization develops for U values of up 
to 1.6 eV, while for U = 2.0 eV a small occupation im- 
balance occurs. The orbital polarization increases with 
increasing /7, and for U > 3.0 eV the system is essentially 
fully polarized. The orbital polarization is accompanied 
by the opening of a gap in the DOS, and the system is 
fully insulating for U > 2.0 eV (see Fig. [sJd). The orbital 
polarization is also apparent from the projected DOS, 
indicating nearly exclusive Px/Py character in the occu- 
pied/unoccupied states, respectively. 

Fig. shows the energy gain due to orbital polariza- 
tion for different orbital order patterns. The notation 
(a, b) indicates the orientation of the occupied orbitals 
on the two different O2 units within the simple tetrag- 
onal unit cell, i.e. a/b can indicate occupation of either 
7T*^y or (tt* iTT*) /a/2 orbitals. All results presented so 
far refer to the {x^y)-type ordering, i.e. an anti-ferro- 
orbital ordering of occupied ttJ and tt* orbitals. Con- 
vergence of the various orbital patterns can be achieved 
by varying the initialization of the orbital occupation 
matrix. As expected from Eq. Q, the calculated en- 
ergy gain associated with orbital polarization (Fig. [4^) 
increases linearly with U for U > 3.0 eV, where the sys- 
tem is essentially fully polarized. Fig. ^ shows the en- 
ergies of the four different ordering patterns relative to 
the (x^y) type ordering [20]. The energy differences be- 
tween the different cases are rather independent of U 



(for U > 3.5 eV) and are all of the order of ~10 meV, 
with the diagonally-oriented anti-ferro-orbital ordering 
{x -\- y^x — y) being energetically most favorable. The 
energy scale of ^10 meV suggests a relatively high or- 
dering temperature of ~100 K. We point out that we 
do not consider any lattice distortion accompanying the 
orbital order so that these energy differences represent a 
purely electronic effect. 

In order to gain further insight into the mechanism 
underlying the orbital ordering, we construct a minimal 
tight-binding (TB) model for the antibonding tt* MOs, 
similar to the one discussed in Ref. [12] for KO2. In 
addition to the kinetic hopping term, we also include 
a mean-field Hubbard interaction resembling the 
correction term of the GGA+/7 approach (Eq. ([T])). The 
hopping amplitudes are obtained by constructing maxi- 
mally localized Wannier functions [21^ 22^ (MLWF) rep- 
resenting the antibonding ttJ^ states around the Fermi 
level. Hopping up to 4th nearest neighbors is considered 
in the TB model (indicated by thick solid/green lines in 
Fig. [1]), leading to excellent agreement between the TB 
and GGA bands (see Fig. More details about the 

TB+/7 model will be presented elsewhere. 

Within this TB+/7 model, we always obtain a polar- 
ization of the (x ± X ± ?/)-type, i.e. an orientation of 
the occupied tt* orbital along the diagonal in-plane di- 
rection. The corresponding ferro- and anti-ferro-orbital 
configurations are degenerate for the AFM case. Fig. [sJd 
shows the evolution of the band gap as function of the 
Hubbard U for both the TB and the GGA+/7 calcula- 
tions. Note that, as previously discussed, a full orbital 
polarization (An ^ 1) of the tt* MOs corresponds to an 
occupation difference An ~ 0.5 in the basis of atomic 
p orbitals used for the GGA-\-U calculations. Thus, the 
same on-site energy splitting Ae = U An corresponds to 
twice the value of U in the GGA+/7 calculation compared 
to the TB model. Within the TB model a fully insulat- 
ing state is obtained for /7 > 0.5 eV, where the width of 
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FIG. 6: (color online). Left: Point charge model for the elec- 
trostatic interaction. Occupied orbitals at each unit are 
approximated by four — e/4 charges located at the positions of 
the MLWF extrema {d = 0.68 A, e = 0.42 A) and +e charges 
are located at the Rb^ sites. The orientations of the two 
occupied MOs are represented by angles a and p. Right: 
Corresponding electrostatic energy per unit cell as function 
of a and p. 



the band-gap follows a linear relation = U — Uq 

with Uo = 0.34 eV. The corresponding expectation for 
the GGA+/7 calculation is thus AE^^^ = U/2 -Uq. It 
can be seen from Fig. [sJd that this relation is well ful- 
filled for U > 3.5 eV, where the system is indeed fully 
polarized (see Fig. The strong suppression of the 

orbital polarization for U < 3.0 eV in the GGA+/7 case 
indicates the importance of electrostatic effects which are 
not included in the TB+/7 model. 

Full orbital polarization results in a pronounced asym- 
metry of the valence charge density, which is unfavorable 
from an electrostatic point of view. In order to assess 
the influence of this effect on the orbital order, we calcu- 
late the corresponding electrostatic energy contribution 
within a simple point charge model (see Fig. [6| . It can 
be seen that the corresponding energy has pronounced 
minima at a, /3 = 45° or 135°, corresponding to occupied 
(ttJ ±7r*)/V^ orbitals. The electrostatic interaction thus 
favors the same diagonal orientation of the occupied or- 
bitals as the hoppings, and furthermore the degeneracy 
between the corresponding ferro- and anti-ferro-orbital 
orientation is lifted in favor of the {x -\- y^x — y)-type 
anti-ferro-orbital ordering (a = 45°, P = 135°). This is 
also consistent with the result of the GGA-\-U calcula- 
tions (Fig.[4]3), even though the very simple point charge 
model is not able to fully account for the energetics be- 
tween all the different configurations. 

In summary, our calculations show that the p elec- 
tron magnet Rb02 exhibits a very strong tendency to- 
wards orbital polarization, driven by strong on-site in- 
teractions. An orbitally ordered insulating state appears 
for rather small values of the Hubbard U (/7o = 0.34 eV 
for the MO basis). We note that U ^ 3.55 eV was de- 
termined in Ref. [12j for the molecular tt* orbitals in 
KO2 based on constrained LDA calculations. The cal- 
culated energy differences between different orbital order 
patterns are of the order of 10 meV, and further analysis 
indicates that the orbital order pattern is determined by 



both hybridization effects (hopping) and electrostatic in- 
teractions. While our study is limited by the rather small 
unit cell size used to explore different orbital order pat- 
terns, our results give clear evidence for the importance 
of correlation effects in p electron magnets. An interest- 
ing question that arises is how the correlation-induced 
orbital order discussed here affects the structural distor- 
tions via Jahn- Teller-like effects and elastic interactions. 
The same effects are likely to be important also for defect- 
related p electron magnetism, which means that simple 
LDA/GGA approaches cannot be used for reliable predic- 
tions of ordering temperatures and other characteristics 
in such systems. 
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